FIRST-ORDER EXPANSION FOR THE DIRICHLET EIGENVALUES 
OF AN ELLIPTIC SYSTEM WITH OSCILLATING COEFFICIENTS 



CHRISTOPHE PRANGE 



Abstract. This paper is concerned with the homogenization of the Dirichlet eigenvalue 
problem, posed in a bounded domain f2 C R 2 , for a vectorial elliptic operator — V-v4 E (-)V 
with e-periodic coefficients. We analyse the asymptotics of the eigenvalues A e,fe when 
e — > 0, the mode k being fixed. A first-order asymptotic expansion is proven for A e ' fc 
in the case when il is either a smooth uniformly convex domain, or a convex polygonal 
domain with sides of slopes satisfying a small divisors assumption. Our results extend 
those of Moskow and Vogelius in |17| restricted to scalar operators and convex polygonal 
domains with sides of rational slopes. We take advantage of the recent progress due 
to Gerard- Varet and Masmoudi |1 H I10| in the homogenization of boundary layer type 
systems. 



1. Introduction 

This paper is devoted to the homogenization of the Dirichlet eigenvalue problem 

, , f -V- A('f)X7v £ = X £ v £ , 

[ ' \ v £ = 0, x G dO, 

posed in a planar domain f2 with periodic microstructure. Some reasons for the study of 
the asymptotical behaviour of the eigenvalues when the period e — > are expounded in 
[19j. Among physical motivations is the analysis of low frequency vibrations in periodic 
composite media. Significant progress in the direction of a better understanding of the 
asymptotics of A e when e — > has been achieved first by Santosa and Vogelius in |19| 
then by Moskow and Vogelius in |17] under weaker assumptions. Our work extends the 
results of |17] to the case of elliptic systems and more general domains £1. Moreover, error 
estimates have been improved. 

Before entering into more details, let us state our mathematical framework. Let N G N, 
N > 1. Throughout this paper, $7 stands for a bounded open subset of M 2 , v £ = v £ {x) G 
1^ and A = A afS (y) G M N (R) is a family of periodic functions of y G T 2 indexed by 
1 < a, (3 < 2. Therefore, taking advantage of Einstein's convention for summation: 

(v.A(f)v^)=4.(^f(f)^). 

All along these lines, C > denotes an arbitrary constant independant of e. The main 
assumptions on A are: 

(Al) ellipticity: there exists A > such that for all f = (£\? 2 ) G 1^ x R N , for all 

y s m 2 , 

A£ Q • £° < A a P(y)£ a ■ f < X- 1 ^ ■ £ a ; 
(A2) periodicity: for all y G M 2 , for all h G I?, 

A{y + h) = A{y); 
(A3) regularity: A is supposed to belong to C°°(J7); 
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(A4) symmetry: for all 1 < a,0 < 2, for all 1 < i,j < N, A"f = A?". 

Unless otherwise specified, we always assume (Al), . . . , (A4). In a very classical fashion, 
boundedness of fi and ellipticity of A imply, through Poincare inequality and Lax Milgram 
lemma, that the linear mapping 

T £ : f £ L 2 (Q)^u £ e H^tt), 

where u £ is the unique weak solution of (1.1) with r.h.s. equal to /, is well defined, 
continuous and injective. If one composes T £ with the compact injection of Hq(Q) in 
L 2 (f2), one gets a compact operator, again denoted by T e , from L 2 (Q) in itself. Assumption 
(A4) tells that T £ is self-adjoint. 

From the previous considerations, we know that our eigenvalue problem (1.1) is well 

posed. There exists a sequence of eigenvalues < A £ '° < A £>1 < . . . X £,k fc ~ > °°> oo and a 
hilbertian basis (v £,k ) of L 2 (Q) of corresponding eigenvectors. To tackle the issue of the 
asymptotical behaviour of (\ £ ,v £ ) we deeply use the periodic structure of the problem at 
microscale contained in (A2). We expand, at least formally, v £ and A e in powers of e 

(1.2) v £ (x)^v°(x, X ) +ev 1 (x, X ) +e 2 v 2 (x, X )+... 



(1.3) X s ^ X° + eX 1 + e 2 \ 2 + ... 

where for all i G N, v l = v l (x,y) is periodic in the y £ T 2 variable. Plugging (1.2) and 
(1.3) in (1.1) and identifying the powers of e yields that v° does not depend on y and that 
(A u , v°) solves the homogenized eigenvalue problem 

, . (-V-A°Vv° = X°v°, xeft 

{ ' \ v°= 0, xedQ ■ 

As usual, the constant homogenized tensor ^4° = A°' al3 € Mn(M) in ( |1.4[ ) is given by 

A o,«0. = f A^(y)dy+ f A^(y)d y y(y)dy, 
where the family x = X^iv) ^ Af?v"(]R), y G T 2 , solves the cell problem 
(1.5) - V y ■ Aiy^yx 1 = d ya A a \ y G T 2 and / X 1 {y)dy = 0. 

Note that ^4° fulfils assumptions (Al) and (A4), so there exists a sequence of eigenvalues 

< A - < A ' 1 < . . . A°> fc — — > oo and a hilbertian basis (v°' k ) of L 2 (Q) of corresponding 
eigenvectors. Let T° denote the operator similar to T £ with ^4° in place of A £ := A Q). 

1.1. What is at stake? Let us now focus on the convergence properties of the eigenvalues 
X £,k when e — > oo and the mode k is fixed. The first thing we know from pi), among other 
papers, is that for all k G N, 

1 1 



(1.6) 
and that 



\e,k \0,k 



^ 1 1 rpe t^O 1 1 

- IK 1 \\C(,L 2 (Q)) 



in £(L 2 (17)) norm. Therefore, for all k G N, 

X e,k e~r^ A L 

no matter wether X°' k is simple or not. When N = 1, i.e. the system (1.1) is a scalar 
equation, on condition that one has enough regularity on v 0,k (y Q,k G H 2 {Q) is sufficient) 
so that an estimate like 

(1.7) ||^(x)-^(x)|| i2(n) <C,||^|| H2(n) 
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holds, one has in addition the error estimate 

(1.8) |A £ ' fc -A°' fc | < C k e. 



Estimate ( 1.8 ) is the starting point of the work of Moskow and Vogelius. Indeed, it leads 



A 



e.k 



X 0,k 



to natural questions, such that: 

(1) What are the limit points of 

(1.9) 

t 

(2) Is there possibly one unique limit point? 

(3) What is the next term in the asymptotic expansion? 

When it does not lead to any confusion, we shall now omit the exponent k: 



A £ 



A 



e.k 



..£,k 



(resp. A 
(resp. v° : 



v 



o,k 
0,k 



A U ' K ) 



In |17j . Moskow and Vogelius get an asymptotic formula for the eigenvalue A e of the 
scalar equation (N = 1), valid up to the order 1 in e, provided that v° G i/ 2 (f2), which 
is actually true for sufficiently smooth domains Q, (convex or C 2 domains). They fully 
describe the first-order corrections, i.e. the limit points of (1.9), in the case when f2 is a 



convex polygonal domain with sides of rational slopes. In this case there is a continuum 



of accumulation points for (1.9). 



Theorem 1.1 (Moskow and Vogelius in |17|). Assume that Q is a convex polygonal domain 
with sides of rational slopes and that N = 1. Assume furthermore that X° is a simple 
eigenvalue. 

Then, there exists #?, E L 2 (£2) solving an explicit homogenized elliptic boundary value 
problem (see (3.6)), and a sequence (e n ) tending to such that 



(1.10) 



A £ « = A + e n X° I r bl (x)v°(x)dx + o{e n ). 
Jn 



1.2. Difficulties and strategy. One faces essentially two kind of difficulties in proving a 



first-order asymptotic expansion for the eigenvalues like (1.10): the first is linked with the 



homogenization of boundary layer type systems, the second has to do with the regularity 
of solutions to elliptic systems in nonsmooth domains like polygons. We sketch how these 
difficulties are addressed by Moskow and Vogelius and how we extend their results to the 
case of elliptic systems and more general polygonal or smooth domains fi. 



1.2.1. Homogenization of boundary layer systems. While v° solves (1.4) with a homoge- 



neous Dirichlet boundary condition on dQ,, v 1 (ji) does not cancel in general on 80,. For 
this reason, the formal asymptotic expansion ( |1.2[ ) is inadequate to establish a first-order 
expansion for the eigenvalues. Boundary layers need to be taken into account. Considering 



v £ (x) w v°(x) + e v l (x, ^+v l b f{x) 



v'lx, ^) +v 2 b f{x) 



+ 



where v % b f solves 



(1.11) 



-V-A(f 



i.e 
v bl 



V. 



bl 



-v \x 



(*.f). 



x G £1 



proves to be more relevant than (1.2). Moreover, ?9f,, appearing in (1.10), comes from the 



homogenization of an elliptic boundary layer system like ( |1.11[ ). 
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The heart of the proof of theorem is subsequently the homogenization of boundary 
layer type systems 

with ip = ip(x,y) := &(y)<po(x), being, unless stated otherwise, a smooth function on T 2 

and ifo G Hz^QQ). 

Compared to u £ solution of 

(1.13) 



v-^(f)Vu £ = /, xen 



u £ = ipo, x 6 dQ, 

whose homogenization is now a classical topic, the analysis of the asymptotics of uf, when 



is complicated by the oscillating boundary data in ( |1.12 ) for at least two reasons: 



(1) We lack uniform a priori estimates for u\i in H l (Q) norm. This is due to the fact 
that 

||^(x,f)|| i, =0( e -§). 

(2) The behaviour of the boundary layer along the boundary d£l deeply depends on 
the interaction between the periodic lattice and the boundary. Thus one can not 
expect in general periodicity of the boundary layer along the boundary. 

The proofs of Moskow and Vogelius intensively rely on a result due to Avellaneda and 
Lin, which addresses a priori estimates for elliptic equations in domains £1 with quite low 
regularity. 

Theorem 1.2 (Avellaneda and Lin in j6] theorem 3). Assume that is Lipschitz and 

satisfies a uniform exterior sphere condition. Assume furthermore that N = 1. 

Then, for all 1 < p < oo, there exists C > such that for all boundary data function 



(/?(•, -) S L p (dQ), there is a solution G L P (Q) of (1.12) satisfying 
( L14 ) KlliP(n) ^ c H-> i) L>(«)) ■ 



Its usefulness for our boundary layer system ( |1.12 ) comes from the following simple 
remark: (/?(•,-) is bounded in L 2 (dQ) norm, but not in H^(dQ) norm. Consequently, 
|| u fj||^2(Q) = 0(1). An estimate similar to (1.14) holds also for elliptic systems, yet under 
stronger regularity assumptions on f2. 

Theorem 1.3 (Avellaneda and Lin in [5] theorem 3). Let N be any positive integer. 

Assume that Q is C ,a with < a < 1. 

Then, the conclusion of theorem remains true. 



This theorem can be applied when f2 is smooth, but due to its strong regularity assump- 
tion on f2, it is useless in the case when Q is a polygonal domain. A precise analysis of the 
boundary layer system is needed in this case. Beyond the results of Avellaneda and Lin, 



the analysis of (1.12) has been carried out in the context of: 

(1) convex polygonal domains f2, first with edges of rational slopes by Moskow and 
Vogelius in [IT], Allaire and Amar in |3j (scalar case), then with edges of slopes 
satisfying a generic small divisors assumption by Gerard- Varet and Masmoudi in 
HT]; 



(2) smooth domains with uniformly convex boundary by Gerard- Varet and Masmoudi 
in the recent paper [10J. 



This recent progress in the homogenization of (1.12), due to Gerard- Varet and Masmoudi, 
opens the way to our generalizations. 
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1.2.2. Regularity. Besides the issue of the homogenization of boundary layer systems comes 
the problem of regularity. Regularity is required in order to carry out the energy estimates 
of the paper. Of course this is only a problem when Q is a polygonal domain; if fi is 
smooth, all functions we deal with belong to C°°(0). 

Assume now that f2 is a convex polygon. In the scalar case the results of Grisvard in [12J 
(theorem 3.2.1.2) and |13j (section 2.7), recalled in fT7], yield that v° £ H 2 (n), because 
of convexity. We even know better. Indeed v° solves (1.4) with r.h.s. X°v° £ if 1 (il). 
Therefore, v° £ H 2+W (Q), with < u. 

This H 2+U (n) regularity on v° appears to be the minimal regularity one has to assume 
in order to get a first-order expansion like (1.10). It is a corollary of the work of Dauge 
[8] on the one hand, and Kozlov, Maz'ya and Rossmann [14J on the other hand, that the 
results of Grisvard extend to systems with constant coefficients. More precisely: 



Theorem 1.4. Let N > 1 and u° £ Hq(£1) be the unique solution of (1.4) with r.h.s. 
equal to f £ H^ 1 ^). Assume that £1 is a convex polygonal domain. 

(1) If f £ H- 1+U1 (n) with < u < 1 and u / \, then u° £ H 1+U >(fi). 

(2) If f £ L 2 (Q), then u° £ H 2 (Q). 

(3) If f £ H 1 ^), then there exists < to < 1 such that u° £ H 2+U1 (£l). 

Let us give a sketch of how to deduce such regularity statements from [8] and |14j (see 
these references for more details). We know from |8] (see lemma 5.7) that for all s > 0, for 
all vertex x £ dO, 

{0 < Re(A) < s} n SpC x is finite, 
where SpC x denotes the spectrum of a pencil associated to our problem at vertex x. 
Besides, Kozlov, Maz'ya and Rossmann prove in |14| . theorem 8.6.2, that for strongly 
elliptic systems, with constant coefficients, satisfying the symmetry assumption (A4), 
posed in the convex polygonal domain 0, 

{0 < Re (A) < l}nSp£ z = 

for all vertex x. This fact collapses if Q has at least one angle > ir. Yet f2 being polygonal 
and convex, it follows from [8], in particular paragraph 7.16, corollary 5.16 and theorem 
5.5, that the operator 

: u £ H s+1 (n) n flj(n) ■— > -v • aQvu g h'-^q) 

is a Fredholm operator for all < s ^ | satisfying {Re(A) = s} n SpC x = for all vertex 
x £ dfl. At this point, one deduces that is a Fredholm operator for all < s < 1, 
s i. Moreover, there exists < ui < 1 such that L^ 1+ ^ is a Fredholm operator. 

Let < s such that is a Fredholm operator and let / G iiP -1 (Q). Two situations 
are possible. If there is a vertex x £ d£l and A G {0 < Re(A) < s} n Sp£a;, then theorem 
5.11 in [8] yields the existence of u? reg £ H 1+S (Q), the regular part, and u° sing £ ff 1+7 (0), 
the singular part, with < 7 < min A6{0<Re(A)<s}n |j SpCx Re(A), such that 



_ .,0 



On the contrary, if for all vertex x, {0 < Re(A) < s} n Sp£ z = 0, then u = nr reg is in 



H 1+S (Q). The two first points of theorem 1.4 as well as the third now easily follow from 
the preceding results. 

Each point of theorem 1.4 plays a role in our reasoning. One can alternatively invoke 
weaker regularity results such as: 

Theorem 1.5 (Agranovich in [2] theorem 1). Assume that £1 is a Lipschitz domain. Let 
u £ £ Hq(Q) be the unique variational solution of (1.1) with r.h.s. equal to f £ 1 (Q) . 



Assume furthermore that f € H 1+t "(fi) with < ui < \. 
Then, u e G H 1+U >(n). 

1.3. Outline of our results. This article answers relevant questions asked by Moskow 
and Vogelius. Quoting [T7] (section 5): 

It should be extremely interesting to derive a similar representation formula 
for polygons with sides of irrational slopes or for smooth domains. In 
particular, it would be interesting to see if this leads to a single first-order 
correction for any eigenvalue. 

We manage to free ourselves from the assumption N = 1. Our main results then sum 
up in the upcoming theorems. We treat separately two different classes of domains Q: on 
the one hand very smooth domains, on the other hand convex polygonal domains. All the 
definitions we use are made rigourous later in the paper. 

Assume that A is an eigenvalue of order m. Let A = A 0,fc = A 0,/c+1 = . . . = A°' fc+m_1 be 
the eigenvalues repeated with multiplicity. We call E\o the finite-dimensional eigenspace 
associated to the eigenvalue A . Note that the eigenvectors v 0,k , . . . , u °> fc + m_1 form an 
orthogonal basis of E^o. 

Our first theorem is concerned with smooth domains f2. 

Theorem 1.6. Assume that Q is a smooth C°° bounded domain with uniformly convex 
boundary. 

Then, for every < j < m — 1, there exists a unique -d* bl E L 2 (f2) such that for all 
< 7 < h 

-i 



(1.15) 



m— 1 



111 f 



1 



m * — ' \ £ ' k +3 
3=0 



\0 m -! r 

X ° + £ —Yl / r i,bl( x ) ■v 0Mj (x)dx + O(e 1+ ^). 
m j=0 J Q 



The next theorem faces the same problem for convex polygonal domains Q. 

Theorem 1.7. Assume that 0, is a convex polygonal domain with sides of slopes satisfying 
a generic small divisors assumption; see section\3A 



(1) Then for every < j < m — 1, there exists a unique •§*■ bl S L 2 (Q) and < 7 such 
that 



(1.16) 



1 



m— 1 



1 



3=0 



X° + e 



m—l „ 

E / ®hi^)-^ k+ i{x)dx + 0{e x ^). 



(2) // E)p C H 3 (£l) n C 2 (J7), then for every < j < m — 1, there exists a unique 



'd* bl G L 2 (Q) such that 



(1.17) 



in— 1 



1 E 

m — s 



3=0 



-1 



\0 m -! (■ „ 

E / r i,bl( x ) ■v 0Mj (x)dx + O(e2). 



We stress that theorem |1.7| has two parts. The first point is a general result: due to 
the assumptions on A (in particular (Al) and (A4)) and on ft (polygonal and convex), 
the H 2+U1 (n) regularity, with < ui, needed on the eigenvectors for the proof, happens 
to be automatically fulfilled. The second part of the theorem states an optimal result in 
view of our proof, in terms of convergence rate, but needs to assume more regularity on 
the eigenfunctions. 

There is a analog of theorem |1.7| in the case of a convex polygon with sides of rational 
slopes, which improves theorem 1.1 Estimate ( 1.16| ) (resp. ( 1.17| )) still holds however up 
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to the extraction of a subsequence (e n ). Throughout the paper, we indicate how to adapt 
the proofs to this case. 



When A is simple, (1.15), (1.16) and (1.17) yield the first-order expansion 



A° + eX 



r bl (x)v°{x)dx + o(e 1+ "'). 



valid for appropriate exponents 7. 

A consequence of these two theorems 1.6 and 1.7 is that the first-order correction to the 
eigenvalue A is identified and unique. Furthermore, it appears in the course of the proof of 
theorem |1.7| that •&* bl is a solution of an homogenized elliptic boundary value problem (see 



(3.7)), whose data can be made explicit. It thus opens the door to numerical computations. 



1.4. Organization of the paper. In section [2] we prove some corrector results for u £ 
solution of (1.13). Such estimates do exist in the litterature, but we focus on minimal 
regularity assumptions. In particular, we extend the bounds of |17j to elliptic systems 
(non necessarily symmetric) and get new ones, which are useful in the rest of the paper. 
Section [3] is devoted to the homogenization of boundary layer type systems. We analyse 
the convergence in L 2 (f2) of "& £ solution of (1.12) with ip(x,y) := — X a (y)dx a v°(x), in the 
two different settings: f2 is a smooth domain with uniformly convex boundary or a convex 
polygonal domain with the additional diophantine condition on the slopes. This work is 
the central step in the proof of theorems |1.6| and |1.7| The final step is done in section [4j 
where a first-order correction formula for A e , in terms of the limit of {} e v bl when e — > 0, is 
obtained. 



2. Some error estimates 
Let / G and ip G Hl(dQ). The solution u e of 



(2.1) 



-V-A(^)Vu £ = f, x£tt 
u £ = (fo, x G d£l 



exists, is unique in ff 1 (Q) and converges strongly in L 2 (S7) towards u° G H 1 ^) solving 
the elliptic system 



(2.2) 



-V • yl°Vn ! 



/, x G 
(fo, x G d£l 



We focus here on estimates in norm showing how fast this convergence takes place. We do 
not need assumption (A4), i.e. the symmetry of A. 

2.1. Multiscale expansions. Before coming to the estimates, let us recall some basic 
facts about multiscale expansions. In the same fashion as v £ (see (|1.2[)), we expand u £ : 



(2.3) 



+ eu 1 1 x 



- e 2 u 2 (x, ^ 



+ ... 



Plugging (2.3) in (2.1) and identifying the powers of e yields, at least formally, 



(1) that u° solves the homogenized system 



-V • A°Vu 



f, x G 
ipo, x G dQ 



(2) that u 1 = u l (x,y) := x a (y)d Xa u°(x) + v}{x) where x° is the function defined in 
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(3) and that u 2 = u 2 (x,y) := T a P (y)d Xa d Xp u° (x) + X a {y)9 Xa u 1 (x) + u 2 (x) where r a/3 
solves 

-V y ■ A(y)V y V a P = B a P- I B^(y)dy, y E T 2 and / X \y)dy = 

Jt 2 Jt 2 



with 



B a f> := A af> + A a W y y + 8 Vl {A^x P ) 



We always assume that u 1 = u 2 = 0. 

The first-order correction u 1 (•, j) to u e does not satisfy homogeneous Dirichlet boundary 

conditions on d£l. It is therefore responsible for a O term in the H 1 ^) estimates 

involving u . In order to correct this, one introduces a boundary layer function $Zm 



solution of (1.12) with ip(x,y) := — u 1 (x,y) = — X a {y)dx a 'U (x). Note that u 1 (•, -) +1?^^ 
belongs to Hq(Q). 

2.2. Error estimates. We extend here the estimates of Moskow and Vogelius (cf. |17j 
section 2) to systems. 

Proposition 2.1. Assume that u° E H 2 (Q). 
Then 

(2.4) \\u £ (x) - u°(x) - eu 1 [x, f ) - ^(^H^ < Ce \\u°\\ H2(n) 

for C > independent of e and u°. 

The proof relies on energy estimates on the error 

e £ := u £ (x) — u°(x) — eu 1 (x, -) — e"D £ u bl {x). 

It is a solution of the following system 

(2.5) 

where 



-V-yl(f)Ve e = r e , xen 
e 6 = 0, x E dQ 



(2.6) 



r e := / + V • 



+ eV- 



A(f)v«'(^f 



We intend to prove that the ff^Q) norm of e e is of order e by showing that the source 
term r £ in (2.5) is of order e in i? _1 (0). It is not clear, looking at (2.6), that the latter 



is true. To face this issue, we invoke the classic key lemma, which can be proven using 
Fourier series expansions: 



Lemma 2.2. Let v 

Assume 



oopm2. to 2 > 



E C°°(T 



V • v = and I v = 
T/ien there exists tp = ip(y) E M suc/i f/iai t> = V _L V = 



5iV ' 



Proof of proposition 2. 1 Expanding the source term r £ yields 
(2.7) 



V y ■ A(y)V y u l 



x, - ] H- 



V y ■ A(y)V x u c 



+ f+ V x -A(y)V x u° 



+ e 



V x • A{y)V x u l 



x 



X. 



V x ■ A(t/)V/j (x, |) + [Vj, • A(y)V :c n 1 



x 



The leading idea is to get rid of terms of order and —1 in e. We call 

v := A(y)V y u 1 + A(y)V x u° 

and notice that 

V y ■ v = V y ■ A(y)V x u° + V y • A(y)V y u\x, y) 

= d ya A^{y)d x ,u° + dy a {A^{y)d yfj x l (y))d x y 
(2.8) = 



because x 7 solves (1.5). Thus the e 1 order term in (2.7) cancels and it remains to handle 
the zeroth order term: 

(2.9) / + W x • A(y)V x u° 



x, — ) + 

£ 



Vx • V 



x,- ) + 

x\ 
*,-) + 



Vj, • A^V^ 1 



x, 



V y ■ (v - A°Vu°) = 



f (v - A°Vu°) = 0. 



Here again, we take advantage of (|2.8l) and the definition of A°: on the one hand 



and on the other hand 



It follows that one can apply lemma |2.2| component by component, and get a function 

ifi = ip(x, y) 6 ~K N such that 

v - A°Vu° = V^. 

Due to the fact that v — A°Vu° is a function of separated variables x and y, ip itself is and 
factors into 

(2.10) ip(x,y) = ^(y)Vu°(x). 

The function ^ = (& a (y)) 1<a<2 ^ M]y(R) 2 is given by the lemma and is therefore of class 
C°°. As u is assumed to be in .£f 2 (Sl), tp is in iif 1 (£7) with respect to x. We set 

w := V^iJj 

of regularity L 2 (J7) towards x, we compute 

VyW = Vy V^ = -V x ■ V^V = -V a • V - / 



and use this equality to simplify (2.9) 
(2.11) / + 



• v 



X' 



.r 



• w 



+ 



V y • A{y)V x v} 



x, - I . 

e. 



Finally, using V x • w = one obtains 



Vy ■ W 



Vy ■ W 



Af-lV-uMx 



x, 













e 




(x,f) +£ V- 


(*>!) 











(2.12) 



-eV 



w\ x 



Af-JV^ 1 (x, 



It remains to estimate the H~ 1 (Q) norm of r £ . The expression (2.12) is convenient for 
two reasons: it is a sum of two terms of order 1 in e and is written in divergence form. 
Moreover, w(-, -) as well as A(-) V^u 1 (•, -) belong to L 2 (Q) and we have 

,011 



||A(-)V :E it 1 (-, -)|| L 2 (Q) < ClKHi^n) 
Consequently, for all £ ^(fi), by Cauchy-Schwarz inequality 

>(., f),«.)>_ 



-eV 



(n),a3(n) 

£ 



10 X 



eV 



e / -wj^x, —J • V(j)(x)dx — e / A^-^jVxU 1 (^x,-^J ■ V(j)(x)dx 



in 

< Ce \\u Q \\ H2(sl) Ikllifi(n) 

which concludes the proof. 

Corollary 2.3. Assume that vP £ H 2 (Q). 
Then 



□ 



(2.13) 



\u e (x) -u°(x)\\ L2{n) 



< Ce* \\u 



oi 



Proof. By the triangular inequality, we get 
\\u £ (x) -u°(x)\\ L2(n) 

< \\u £ (x) - u°{x) - eu 1 (x, f ) - e$ £ uM {x)\\ Hl{n) + e\\u 1 (x, f ) || L2(n) +e ||^W|| L2(n) • 

The estimate (2.13) now follows from ( |2.4[ ), the uniform boundedness of ti 1 (-, -) in L 2 (r2) 
and from the bound 

( 2 - 14 ) ¥U\» {Q) < ||^||fli (0) < C ll^(f)^^ Q (^)|| H | (an) < Ce-\ V\\ H2{n) . □ 



From corollary 2.3, one easily gets a similar L 2 (Q) estimate under weaker assumptions 



on u°. 



Corollary 2.4. Assume that u° G H 1+LU (tt), with < oj < 1. 
T/ien 



(2.15) ||u £ (x) - u°(x)|| L2(Q) - wo- (r ii^i+c^) 



II h_n (it j 

Proof. A straightforward energy estimate on the elliptic system satisfied by u £ — vP yields 
(2.16) \\u £ (x) - u°(x)\\ L2{n) < C \\u £ (x) - u (x)\\ H , {n) < C Ik 01 
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Inequality ( 2.15[ ) now comes from interpolating (2.16) and (2.13). The main idea is to 
think of ( 2,16[ ) (resp. (2.13)) as the statement that the linear operator 

u° .— >• u £ {x) - u°(x) 

is bounded from i^ 1 (SI) to L 2 (f2) (resp. from H 2 (Q) to L 2 (S7)) and then to interpolate. □ 

We call now the solution of (1.12) with tp(x,y) := —u 2 {x,y) 1 whose introduction is 
motivated by the same reasons as "Q £ u and state the proposition: 

Proposition 2.5. Assume that u° G H 3 (Q). 
Then 



(2.17) 



: (x) - u°(x) - eu 1 (x, j) - ei9 £ uM (x)\\ L2{n) < Cel ||«°|| HS(n) 



Proof. The proof of (2.17) relies on a global energy estimate found in section 3.3: 
(2.18) u £ {x)-u°(x)-eu 1 (x,f) - £ tf s uU (x) - e 2 u 2 (x,£) - s 2 $ 2 / bl (x) 



H 1 (Q) 



0(e 2 ) 



It requires u° G H S (Q,) and it can be showned using the same ideas than those involved in 
(2.4), the key being again lemma 2.2 Following the lines of |llj it becomes clear that the 
precised estimate 



u £ (x) - u°(x) - eu l (x, f ) - e$l bl (x) - e 2 u 2 (x, §) - e 2 ^ £ bl {x) 



holds. It is nothing but a consequence of the way the involved functions factor in the 
product of a function depending only on y and of Vit° (cf. ( |2.10 )). We conclude, applying 
the triangular inequality, that 



\u £ (x) — u°(x) — eu 1 (x, f ) — ed e u w (x)|| 



L 2 {n) 



< 



u £ (x) - u°(x) - eu^x, f) - e$ £ uM {x) - e 2 u 2 (x, §) - e 2 ^ £ bl (x) 

2 



+ e 2 \\u 2 (x, 



e 



u.bl 



# 2 :Ux 



L 2 {n) 



The estimate (2.17) now follows from (2.18), the uniform boundedness of u 2 (-, -) in L 2 (Q.) 
and the (2.14)-like bound 



2,e 
u.bl 



< Ce 2 . 



□ 



We conclude this section focusing on a (2.17)-like estimate for u° satisfying a weaker 
assumption. 

Theorem 2.6. Assume u° G H 2+W (£l), with < u < 1. 
Then 



(2.19) 



\ u ^ x) - u Q {x)-eu\x^) -e^ M (x)\\ < Ce 1+ ^ \\u°\\ H2+u{n) 



As in the proof of corollary 2.4 estimate ( |2.4[ ) (resp. ( 2,17| )) states that the linear 
operator 

u° i — > u £ (x) — v°(x) — eu 1 (x, — ) — bl {x) 



is bounded from H 2 (Q) to L 2 (f2) (resp. from H S (Q) to L 2 (Q)). By interpolating between 
the two linear operators, one gets ( |2.19 |. All details can be found in |17| and apply without 
any change to the case N > 1. 

n 



3. HOMOGENIZATION OF BOUNDARY LAYER TYPE SYSTEMS 



Throughout this section we are interested in the homogenization of the boundary layer 
type system 

r , n / -V-A(f)W^= 0, xesi 

[ ' I K,u= -x a {l)d Xa u°(x), xedn 

that is in the study of the asymptotic behaviour of the sequence bl when e tends to 0. 
This means we both look for a possible limit of the sequence and for estimates in norm 



of the speed of convergence. For all this section, we assume that u° solves (2.2), with 
/ G L 2 (n) and ip = 0. 



This is a crucial step in the proof of theorems 1.6 and 1.7 As explained in the introduc- 
tion, there is no regularity issue when ft is smooth. On the contrary, when f2 is a polygon, 
we concentrate on minimal regularity. That is why we give two convergence rates: the 
first under the minimal assumption u° G H 2+U} (Q,) with < U) < 1, the second under the 
stronger regularity assumption vP G H 3 (£2) n C 2 (fl), where we focus on improving the 
speed of convergence. 

For notational convenience, let us write in this section #f, instead of ~& £ ubl . 

3.1. Smooth uniformly convex domains. Assume that n C I 2 is a smooth (say C°°) 
uniformly convex domain i.e. all principal curvatures are bounded from below; see |7] 
section III. 7 for another definition. The regularizing properties of elliptic operators in 
smooth domains yield that u° G C°°(Q) (see |T] theorem 10.5). Therefore, the boundary 
data function tp(x,y) = —u (x,y) = —X?*(y)d Xa u(x) is a smooth function. Note that we 
do not need assumption (A4), i.e. the symmetry of A. 

Theorem 3.1 ( Gerard- Varet and Masmoudi in [E]). For all 1 < p < oo there exists 
ip* £ LP{dVt) such that •&% converges in L (fi) towards di, G L P {£1) solution of 

-V • A°V0J, = 0, xen 



$* bl = tp*{x), x G dn ' 
Moreover, for all < 7 < jj, 

(3-2) K-*«ll^ (n ) = 0(^). 

We do not attempt to weaken the regularity assumption on f2, which itself implies strong 
regularity on u°. For details concerning the proof and relevant remarks, we refer to [10] . 



3.2. Convex polygonal domains. Let us assume f2 to be a bounded convex polygonal 
domain with M edges, supported by the lines K k of unitary inward normal n' £ S 1 , Thus 



M 



Q = C\{x, n k ■ x > c k ) 



k=l 

with c k G R, and for all 1 < k < M, 

K k = {x, n k -x = c k ). 
Beyond this first assumption on Q we require either 
(RAT) rationality: for all 1 < k < M, 

(3.3) n k G IRQ 2 

or 
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(DIV) small divisors: there exists C, I > such that for all 1 < k < M, 
(3.4) V£ = (6,6) G \ {0}, |P fc x(OI > Ciel"' 



where P fc ± is the orthogonal projector on n* 



As f2 C M 2 , condition (3.4) boils down to 

(3.5) V£gZ 2 \{0}, |n fc -e|>C|cr' 

where • £ := n^i + n k £,2- Note that a vector n € M 2 cannot satisfy both (3.3) and (3.4) 
or fl3~5j ). 

Keeping in mind that i?f; solves (3.1 ), one has the following convergence theorems. Note 
that as soon as we invoke the regularity theorem 1.4 in the proofs, we need the symmetry 
assumption (A4) on A. 

Theorem 3.2. Assume Q satisfies (RAT). Assume furthermore that u° G H 2+U) (Fl) with 
< lj < 1 (resp. u° G H 3 (fl) n C 2 (U)). 

Then there exists a sequence (e n ) and (V k,a '*)i<.k<M G Mtv(M) 2xM such that solution 

l<cK2 

o/ (1.12) converges in L 2 (f2) towards $?, solution of 

o, ie!l 

-y^d^Ofa;), iE90n A" fc , /or all 1 < k < M ' 



(3.6) 



-V • A°V$* bl 



A* 



Moreover, we have the following convergence rates: 

(1) ifu° G H 2+W (Q), then there exists < 7 < | suc/i t/iai 



I ?9 e " 



6/llL 2 (f7) 



(2) G # 3 (ft) n C 2 (fi) 7 then 



0«); 



0(el). 



Theorem 3.3. Assume Q satisfies (DIV). Assume furthermore that vP G H 2+U) {yt) with 
< lj < 1 ("resp. u° G # 3 (ft) n C 2 (U)). 



Then there exists (V k,a, *)i<.k<M G Mn(M.) 2xM suc/i £/ia£ 1?^ solution of (1.12) converges 



Ka<2 



in L (Ft) towards solution of 



(3.7) 



-V • A°V$* bl 



0, 

-V fc ' a '*0 Xa u°(a:), lEffifl if fc , for all\<k<M ' 
Moreover, we have the following convergence rates: 

(1) G H 2+UJ (Q), then there exists < 7 < ~ suc/i t/iai 



(3.8) 



?9 £ 







(2) ifu° G # 3 (ft) n C 2 (n), then 



(3.9) 



bl\\L 2 (n) 



i9* II 



0(ei). 



The only, but major, difference between theorem |3.2| and 3.3 is that, in the small divisors 
case, convergence holds for the whole sequence, whereas in the rational case, convergence 
takes place up to the extraction of a subsequence (e n ), the constant matrices V k)Ci '* de- 
pending on (s n ). 



3.3. Proof of theorem 3.3 The proofs of theorems 3.2 and 3.3 follow the same steps. 
They differ mainly in one intermediate result, which explains why in the rational case, the 
convergence result is true only up to the extraction of a subsequence. Although we focus 
on the (DIV) assumption, we underline the difference with assumption (RAT). 
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3.3.1. Existence of the boundary layer tails. Let us show the existence of the matrices 
V k ' a '*. Let 1 < k < M and 1 < a < 2. We are interested in the boundary layer profile in 
the vicinity of vertex k. Thus, introduce v k j a ' £ solution of 



-V, • A(y)V y v k b i a > £ 



0, 



v 



k,a,e 
bl 



x a (y), yedtl 



fe.s 



where Q k,e : = jy, n k ■ y 



'- > 0}. Let M k E M2QR) be an orthogonal matrix, mapping 



e2 := 



to n k . 



The following theorem describes the profile of V 



k,a,e „, fc > a > £ 



bl 



Theorem 3.4 (Gerard- Varet and Masmoudi in Assume that Q satisfies (DIV). 

Then, 

(1) for all e>0, there exists V k ' a > £ E C°° (j x ] oo \) ; 

(2) there exists a matrix V k ' a '* E M N (R) such that for all f3 E N 2 , for all m E N, there 

k 

is a constant Ci«i m > satisfying for all e > and Z2 > — , 



(3.10) 



c fc 


m\ 




) sup 




7 ziSK 



<c, ; 



\,m- 



Remark 3.5. Note that (3.10) is true not only for m E N but for m E R, m > 0. Indeed, 
[m] denoting the integer part of m, we have 



Z'2 



Z-2 



h ra 

C 

£ 

„k \ m 



< 



Z'2 



H+i 



if 



Z2 



> 1 



£ 



< 1, if 



^2 



< 1. 



3.4 



in terms of w^' a,e instead 



Remark 3.6. If we rewrite the second statement of theorem 

of V k,a ' e we get: for all f3 E N 2 , for all m E N, there is a constant Ci^i m > satisfying for 
all £ > and y E r2 fc,£ , 



(3.11) 



1 + 







yn k 





) dl{v k b r{y)-V k ^) 



< Q 



\P\,m- 



Remark 3.7. If instead of (DIV) we assume (RAT), the boundary layer tails V k ' a ' £ still 
exist. Furthermore, an equivalent of theorem 3.4 states: there exists a sequence (e n ) (here 
lies the main difference between the two assumptions), a constant matrix V k > a '* e Mjv(M) 
such that for all m E N, for all Z2 > — , 



1 + 



22 



) sup (y fc ' a ' e "(zi,z 2 ) -V k ' a '*) 
' Zl eR 



<C\ 



The latter is sufficient to get our results. However, Moskow and Vogelius in [17j, as well 
as Allaire and Amar in [3j manage to prove an improved result under assumption (RAT): 
the convergence of the boundary layer towards its tail is exponential. 

Assume from now on that O satisfies (DIV) . Let < oj < 1 be fixed. The assumptions 
u° E H 2+U (Q) with < to < 1 and u° E H 3 (Q) n C 2 (Q) are treated in parallel. In both 
cases, by Sobolev injection, u° E C 1 (J7). 
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3.3.2. Well-posedness of (3.7). It is enough to prove that the boundary function of (3.7) 
belongs to H^{dVL). One constructs a lifting (f)T, of ipt, := — V k ' a, *d Xa u°(x). There exists 
G = (G 1 , G 2 ) G M N (R) x M N (R) such that 



(3-12) <f>* bl 

Following one can show: 

Proposition 3.8. If vP G H 2 (9), then 
then (jf bl belongs to H 2 (Q). 



G a d T „ 



G Ifu°, in addition, belongs to H 3 (ft), 



3.3.3. Sketch of the proof of the estimates (3.8) and (3.9). Our strategy is to split the 
problem of estimating — into three easier ones. For this purpose, we introduce $fj 
solution of 



(3.13) 



-V-A(f)W 



£,* 

bl 

if'* 
u bl 



0, x G ft 

-V k ' a '*d Xa u°(x), x G 8ft n if fc , for all 1 < jfc < M 



to get via the triangular inequality: 



i9 e 



?9* II < Ili9 £ '* 



? 9* _i_ 



i9 £ '*ll 

V bl llL2(f2) ■ 



Note that #,! is well defined because of proposition 



3.8 



The study of the first term seems to be more classic as the boundary data function 
of (3.7) and (3.13) is not oscillating. The second term, on the contrary, requires a deep 



knowledge about the homogenization of boundary layer systems. 

In fact tff, - $ £ b f is the solution of JTTl2| ) with <p(x,y) = -(x a (y) ~ V k ' a >*)d Xa u°(x), 
for all x G dQ n K k , for all y G M 2 ; we call uj/ e the difference 0£, - It comes from 
that ip defined like this is in H?(d£l). Let v b { £ := - {v k { a ' £ - V k > a >*) d x „u° . 

>M k,£l 



3.8 



proposition 

We expect u b f to be close to 



l,e| 



Il 2 (c) 



< 



M 



fe,e / 



x, 



fc=i 



+ 



^0*0 -EifciiWw 6 ^!) 



M 



L 2 (Q) 



L 2 (H) 

The rest of the proof is thus devoted to estimate each of the terms in the r.h.s. of: 

M 



(3.14) 



i9 e 
\ V U 



?9* II < Ili9 £ '* 



fe=i 



V bi [ X i f ) 



L 2 (n) 



+ 



1.5/ 



(^-E£i«6 e (x,f) 



fc,5 / 



L 2 (a) 



3.3.4. First term in the r.h.s of (3.14). We resort to corollary 2.4 to estimate this term. 
In order to get some convergence rate, we need to have a little more regularity on than 

the lifting (f) bl of the boundary data of (|3.7|) 



3.8 



$ll G H 1 ^). According to proposition 

belongs to H 1+ "(n) (resp. H 2 (n)), provided that u° belongs to H 2+UJ (n) (resp. H 3 (nj). 

Let us treat the two assumptions on u° separately. If u° G H 2+UJ (Q), it follows from the 
first point of theorem 1.4 that has i? 1+7 (il) regularity for all 7 such that < 7 < to 
and 7 7^ g- Therefore, 

^IL 2 (n) = °(^)- 

yields that i9^ G H 2 (Q). Applying corollary 



If u° G -ff 3 (Q), the second point of theorem 
2.3| implies 



1.4 



0; 



M|lL 2 (f2) 
15 



o( £ l). 



3.3.5. Second term in the r.h.s of (3.14). By linearity of the equations, the boundary layer 



tail V k '*(x) of i>w' e (x, •) is equal to 



V k '*{x) = -V k ' a '*d Xa u°(x) + V k,a '*d Xa u°(x) = 0. 



We deduce from theorem |3.4| for all m G N, there is a constant C m > such that for all 
e > 0, for all 



(3.15) 



1 + 



\x ■ n — c \ 



V bi (*> f , 



The uniformity in x comes from the fact u° G C 1 ^) and from the boundedness of £1. 



Proposition 3.9. For all 1 < k < M , 

Proof. Let m G N. From ( ]3.15| ) we get 



k,e i 
Vj [X 



' e) 



L 2 (n) 



o( e l). 



k,e i 



o.r {x, f) 



L 2 (n) 



< c 



jdx 



< c 



where :=*M fc fi 



Therefore, we have to focus on 



duo = £ 



2m 



'[0,oo[ (l + f )" 
For 2m > 1 the integral is convergent and 



[0,oo[ (£ m + V% 



-^dU2 



[o,oo[ (e m + u™ 



We immediately deduce that 



which yields the result. 



jdit2 < 



[0,6] 



T 2m 



+ 



,2m 



k,s 



(*.f) 



L 2 (n) 



[e,oo[ a 2 



0(e) 



du 2 = 0{e- 2m+1 ). 



□ 



Remark 3.10. It is easy to adapt the proof of proposition |3.9| to get: for all 1 < k < M, 
for all 1 < p < oo, 



V bi (*> I) 



LP(H) 



O(ep). For p = oo it is (3.15) with m = 0; for 



1 < p < oo the proof follows the lines of the case p = 2, except that one has to replace 

2 by p. In the same manner, it is very straightforward to deduce from (3.11) that for all 
1 < P < oo, for all G N d , for all m G N, 

i , 



(3.16) 



Q.8 ( v k,a,c(-j \rk.o.i 



V \ u bl 



'(i)-v k 



0(ep). 



3.3.6. Third term in the r.h.s of (3.14). We proceed as usual by carrying out energy 
estimates on the error 

M 

e . _.l,£/'_\ _.&,£ f_ x ^ 



'hi 



U 



bl 



fc=i 



It solves the system 



-V- A(f)Ve; 



r w x ^ n 



£ 

bl 

e li = ¥>w. 
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where the source term 

M 



'bl 



k=l 



+ 



V x ■ A(y)V y v 



X, 



and the piecewise defined boundary function 



PbUdnnK* 



X 



x M 

V k ^)d Xa u\x)-Y,v 
' k'=l 



k'e( x 



M 

E 

k'^k 



k',ef X 



We estimate separately rf, (cf. lemma 3.11) and (pf, (cf. lemma 3.12). 



Lemma 3.11. The source tcvTfi Vfoi is 



(1) of order 0(e?) m H' 1 ^) for all < 7 < f i/ ii° G H 2+U} (Q); 

(2) 0/ order O(el) zn fl-^Jl) z/u° G C 2 (H). 

Proo/. Assume u° G H 2+UJ (Vt) (resp. n° G C 2 (H)). Let 1 < k < M be fixed and consider 



r £ b f : = V - ( .!(- )V,r 



fe.s / 



We focus on the first term of r e h f . For all G -ffg (Q), 



fc, e / a; 



V- A - V x v?f x, 



,(t>{x) 



A 



V x v^f[x, J V4>{x)dx 



< 

< C 



At this point we need to estimate 



L2{n) W\\Hi(n) 
. As 



L 2 {n) 



the idea is to bound the L 2 (ft) norm of this term using a Holder inequality and (3.16). 
Doing so, one has to pay attention to the regularity of u°, and to carefully choose the 
L P {Q) spaces involved. 
If u° G C 2 (U), 



fe,/9,£ 



(l)-v k ^*)d Xa d xpU ° 



L 2 (n) 



< 



k,/3,e 



(x\ _ yk,p,* 



L2(Q) ll^^ n °IL-(Q)- 



The assumption u° G C 2 (£l) plays here the same role as u° G C 1 (il) for (3.15). Use (3.16) 
with p = 2 to conclude. 

If u° G H 2+U) (Q,), we cannot proceed as above because d Xa d X/3 u° does not belong to 
L°°(f2). By the Sobolev injections, is continuously embedded in L q (£l) for all 

1 < q < i^j- Yet d Xa d X pU° is in Take 2<g< r ^ J andg>2 such that 



1 + i = |. Necessarily ^ < Q- Holder's inequality yields 



q q 



L 2 (n) 



Apply now (3.16) with p = q to get 



< 



k,/3,e. 
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L 2 (n) 



o( e t). 



The second term of r^' fc needs to be treated differently. The key ingredient is Hardy's 
inequality: for all <f> G Hq(Q) 



d(x,dn) 



L 2 {fl) 



< \\V4>\\ 



L 2 (Q) 



where d(x,d£l) is the distance from x to dil. Let (ft G Hq(Q). For all 2 < g, q < oo such 
that I + k = 2, 



x, — 1 4>{x)dx 



< Ce 

< Ce 



c k \ 




4>(x) 






L 2 (Q) 


d(x,an) 


L 2 (n) 



U Vp U bl \e) 



L8(H) 



\d Xa d Xrf u°\ 



Li(n) lr ^IIl 2 (^) • 



If u° G H 2+UJ (Q), then take 2 < q < and apply (3.16) with p = q and m 
u° G C 2 (0), then take q = oo and q = 2 and apply (3.16) with p = 2 and m = 1. 

Lemma 3.12. T/ie boundary function (p^ is 

(1) o/ order 0(e w ) in W l ~p' p (dtt) for alll<p<2, if u° G H 2+LU (n); 

(2) o/ order 0(e) in PF 1- *' P (<90) for alll<p<2, if u° G H 3 (tt) n C 2 (n). 



1. If 

□ 



Proof. First of all, using proposition 3.8 one notices that tpL belongs to Hz(dQ). As ip^ 
factors into V( L )Vu° we immediatly get the very rough estimate 



I ^ W H#2 (an) 



o(e-r 



which is far from being enough. We do not try further to get a bound in H 2 (dQ). 

We refer to |H| for the case when u° G H 3 {n) n C 2 {U). If u° G H 2+UJ {VL) the proof 
follows the same scheme, with differences due to the weaker regularity assumption on u . 
Assume for the rest of the proof that u° G H 2+UJ (Q). The edge estimate goes on as in the 
case u° G H 3 (n) n C 2 (n) and one gets for all 1 < p < 2, m G N 



\tjjip £ 



W 1 p' p (dfl) 



0(e n 



where ip is a smooth function on dVL compactly supported in dVLC\K k outside the vertices. 

Let us now focus on the estimate near a vertex O lying at the intersection of K l and K 2 . 
We introduce polar coordinates r = r(x) and 9 = 9(x) centered at O and use a smooth 
function ip on d£l compactly supported in a vicinity of O. Let 1 < p. The tame estimate 



(3.17) 1-1 



w p''{dn) 



<C 11/ 



lL°°(ao) 



3 



+ 5 



L°°(<9!T2) 



l/l 



W p' p (9Q) 



holds for all f,g G L°°(dft) n W 1 "^^). Taking advantage of the fact that H 2+UJ (n) 
injects in C^H), one knows G L°°(dn). Besides, ^ belongs to W 1 '^) for all 
1 < p < 2. Therefore G L°°(dn) n W 1_ p' p (dn) for all 1 < p < 2 and <Kl7\ yields 



II^V&I 



w p' p (dn) 



< 



c{\\^v(,)\\ w ,- h . 



(9n) 



L°°(sn) 



W P' p (dQ) 



18 



By estimating first on dfl D K 1 then on dft fl K 2 one obtains for all 1 < p < 2 

ll^(i)L~ (an) = o(0 

(3.18a) II^Vf^L,,^, = 0(e u 



(3.18b) 



\tpr bJ V{ 



\LP(dn) 



0(e 



->)■ 



Interpolating ( |3.18a ) and ( 3.18b[ ) gives 

< c U r "v(-)\\ i w i {dn) \\^v(-)\\ 



\^V{-)\\ ,.i 



lp (an) 



O e 



u-l+i 



Finally, H^Vw 



0(e w ) which concludes our proof. 



□ 



3.11 



We conclude this section by expounding how to deduce a bound on e| ; from the lemmas 
and 



3.12 We focus on the case G H 2+UJ (n), as the is a little more 

subtle than in the case u G i/ 3 (f2) n C 2 (£7 ), It i s very straightforward to adapt the proof 

in the latter case (see also |11|). Lemma 3.12 gives bounds for <pf, in W 1 p ,p (d£l) for 



1 < P < 2. As we lack an estimate in Hv(dVt), we cannot bound e?, in ff 1 ^). We thus 
have to use results on elliptic equations in divergence form and with source term in some 
VF -1,p (fi) space. Let us state a general theorem that suits to our framework (for references 
see below). 

Theorem 3.13 (Meyers). Let C R d be a Lipschitz domain, A = A af3 {y) G M N (R) a 
family ofC°°(Q) functions. Assume the ellipticity of A. 

There exists a po < 2 such that for all f G -ff _1 (0) if u G Hq(Q) is a weak solution of 
-V • AVu = f in fl-^fi) and if for all p < p < 2, f G W' 1 ^), then u G Wq' p (Q) and 
there exists C(p) > 0, 



W-!.P(f2) " 



Such an estimate originally appeared in the work of Meyers |16J . where the case of 
smooth C 2 domains f2 and scalar equations is treated. It has been extended by Gallouet 
and Monier in [9] to domains f2 with Lipschitz boundary. In their recent survey article 
[15], Maz'ya and Shaposhnikova give very general estimates working for Lipschitz domains 



Q, and systems of elliptic equations. Our theorem 3.13 happens to be a very special case 
of theorems 1 and 2 in [15] . To make the link obvious take m = 1, I = N, a = 0; then 
8 = 1 - ±, W™> a (Q) = W^P(n) and V*(tt) = W 1,p (ft). 

It is important to notice that po (resp. C(p)) only depends on the coercivity constant 
of A (resp. on the coercivity constant of A and p). This makes the theorem applicable to 
our homogenization problem. We know from the proof of lemma 3.12 that 

fbl W 1 '^^). 

l<p<2 

Thus there exists a lifting 0|, of ip^ belonging to W 1,p (tl) for all 1 < p < 2 such that 



\0li\\ w i, P{n) <C(p)\\ip £ bl \ 



w 1 p' p (dn) 



with C(p) independent of e, as usual. The difference e 

'X 



G Hq(Q) solves 



-V -A 



bl 



r £ bl +V-A - Vtf 



As F?, belongs to W 1,p (ft) for 1 < p < 2, we have 



(3.19) 



F 



bi\\w- 1 ^p{n) 



<C(p) 



(dn) 
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Let pq < 2 given by theorem 3.13 Then, for all pq < p < 2, 



-hi 



Let < 7 < |. 
Po < P < 2, 



Then, it follows from (3.19) and from lemmas 



3.11 



and 



3.12 



that for all 



\F £ \\ 



To get an L 2 (J7) estimate on ef, use the Sobolev injection of W l ' p {^l) in L 2 (Q) and, once 
again, our W 1 p' p (dfl) bound on <^| ; . 

4. A FIRST-ORDER ASYMPTOTIC EXPANSION OF THE EIGENVALUES 



This section is concerned with the final step of the proof of theorems 1.6 and |1.7| Let 
E)p be the finite-dimensional eigenspace associated to A . From the ideas explained in the 
introduction, and in particular the third part of theorem 1.4, we know, in any case, that 
E^o C ,ff 2+a; (f]), with < ijj. When Q is a smooth uniformly convex domain, we take 
uj = 1. 

We have recourse to the ideas involved in [17J to prove the asymptotic expansion of the 
eigenvalues. Moskow and Vogelius use abstract estimates due to Osborn in [18]. We recall 
the estimate we need in terms of T £ and T°. Assume that A is an eigenvalue of order m. 
Then, E x o is m-dimensional. Let A = A°' fc = A 0,fc+1 = ... = \°,k+m-i assoc i a ted 
eigenvectors v 0,k , . . . , v °> k + m - 1 form an orthogonal basis of E^o. 

Theorem 4.1 (Osborn in [18]). There exists a constant C > such that 



(4.1) 



1 1 v-im-l 1 1 v- ,m— 1 

\0 ~ m ?-~>j=0 \£- k +3 ~ m ^i=0 



£7=o 1 ((T £ - T°)v°> k +i , 



< C\\{T £ -T° 



where T £ and T° are seen as operators acting in L 2 (S7) and (•, •) denotes the scalar product 
mL 2 (n). 

This theorem is a straightforward corollary of theorem 3.1 in |17| . Its proof really uses all 
properties of the operators T £ and T°, among other things selfadjointness and compactness. 
The first thing to do is to estimate || (T £ — T°)\e x0 || . Let / £ E^o; we call u £ := T £ f and 



u 



T°f. We need to estimate 



u 



u 



\L 2 (n) 



. We can improve the bounds of section 2 



such as ( 2,13[ ). Those bounds are not enough to deduce from (4.1) a first-order asymptotic 
expansion of A £ . The loss of O 



in estimate (2.13) is due to the bad bound of "& £ ub i in 



Ii9 £ II < 
\ u uM\\L2(n) - 



l^wllifMn) - C W xa 



u°(x) 



\H 2 (n) 



If Q, is a smooth domain, (1.7) can be shown thanks to the results of Avellaneda and 



Lin. Theorem 1.3 yields indeed 

ll?9 £ II 
\\ w u,bl\\L 2 (n) 



<C\\<p( 



L 2 ^) - c !I u °IIh 2 (q) • 



The assumption u° € H 2 (£~l) being clearly fulfilled as u° = T°f = -Aj/, it is easy to adapt 
the proof of corollary |2.3| to conclude that: 

(4-2) \\u £ - u°\\ L2{n) < Ce ||« || H 2 (n ) < Ce \\u°\\ H2+0J{n) ■ 

Assume now that is a polygonal domain satisfying either (RAT) or (DIV). A uniform 
bound in e of i9 £ hl does not follow from the results of Avellaneda and Lin. The estimates 
of section [2] are sufficient to get the convergence of the boundary layer in section [3] up 
to t he extraction of a subsequence under assumption (RAT). Actually theorem 3.2 (resp. 



3.3) implies that: there exists a sequence (e n 
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such that 



Lq£n II 

I uM\\L 2 (n) 



\H 2 +"(U) 



(resp. for all < e, ||^m|| l2 (^) < C || u °||^2+^(q))- We conclude, as in the case when Q, is 
smooth, that (4.2) holds. In order to avoid extracting subsequences we now omit the case 
of polygonal domains under assumption (RAT). 

It remains to bound || u °||#2+w(q) by H/H^amy Taking advantage of the equivalence of 
norms on the finite-dimensional space E\o C H 2+UJ (£l) C H 2 (Q), there exists < C such 
that for all w £ E^o, 



(4.3) 



Therefore, combining (4.2) with (4.3), we get 

1 1 T 6 f f" 

which shows that 

\\(Tt-T°)\ Ex0 \\ <Ce. 



Our final goal is to prove (1.15), (1.16) and (1.17). The reasoning, in every case, follows 



the lines of [17] . Estimate (4.1) now sums up in: 

m—l 

I I 

(4.4) 



1 I 

A m 



m—l 



E J^TJ = I E - TV' fe+ V°^) + 0(e 2 ) 



3=0 



3=0 



v 0,k+j jV 0,k+j^ 



Let us focus on i YJj^o (( TS - T°)v°' k+ i , v°' k+ i) and work on <(T e - T°) 
j being fixed in {0,. . . , m — 1}. We ca ll u £ ,k+ i := T £ v°' k+ - This function solves (2.1). 
According to estimate ( |2.19| ) of theorem plij as v°' k+j G H 2+UJ (n), 



where $% M+jM 



(4.5) ( {T £ - T u )v 



solves (3.1) with t)°> fc +J instead of u°. Cauchy-Schwarz inequality implies 

1 



0\„,o,k+j „,o,k+j\ 



A° 



v°' h+j (x) -u £ ' k+1 {x) ) {x)dx 



,e,k+j 



,0,k+3i 



e 

A° 



X a (f) d Xa v > k+ i(x)-v°> k+ i(x)dx + ^ / V £ V!k+jM (x) ■ v°^(x)dx + 0{e 1 ^ 



We intend to show that the term involving \ a m (4.5) is of order 0(e 1+ a). In order to 
carry out integrations by parts, we introduce, for each 1 < a < 2, a periodic C°° solution 
b a = b a (y) E Mjv(M) to 



AX = X a , 



£ 

A" 



the Fredholm property being satisfied as fj 2 X a {y)dy = 0. An integration by part gives 
X a (-)d Xa v > k +i(x)-v°> k+ i(x)dx 

d Xa v°' k+j (x) ■ v°' k+j (x)dx 

' ha (?)) ■ v ( d *y ,k+j (. x ) v °' k+j ( x )) dx 



£ 

A° 



£ 2 A(F /r 



£■ 



eV 



<Ce 2 ||eV(6 a (f))| 

< C£ 2 . 
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We deduce from ( |4.4[ ) and ( |4.5[ ) that 

m— 1 



1 



1 



1 



j=0 



1 m— 1 

- (i T£ ~ T°)v°' k+j ,v°' k+j ) +0{e 



1 



j=0 
m—1 



-T- 

m ^ A 

j=0 



The results of section [3] now apply, in particular theorems 3.1 3.2 (in this case up to 
the extraction of a subsequence) and 3.3 and yield that 



|,Q£ 

\ u v,k+jt>l 



0* 



'v,k+j,bi\\i2(Qj = 0(e y ) 
for suitable exponents < 7: 

(1) for all < 7 < jj, when O is a smooth uniformly convex domain; 

(2) for all < 7 < % (resp. for 7 = h), when is a convex polygon satisfying either 
(RAT) or (DIV) and E x o C H 2+U) (n) (resp. E x o C H 3 (n) n C 2 (H)). 

Therefore, 

m— 1 



1 

A" 



1 m ~ 1 1 

1 v 1 

m ^ A £ ' fc +^ 
3=0 



1 



mA° 

j=0 



K MjM (x)-v ^(x)dx + O( £ L ^). i 



with 7 given above; from the convergence of the eigenvalues \ e < k+ i towards A ' fc+J , we 
deduce that 

-l 

A 



m— 1 



1 1 



m * — ' A e 

j=o 



m 



lit _L p 

3=0 Jn 



which achieves the proof of theorems 1.6 and 1.7 
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